Corotating and irrotational binary black holes in quasi-circular orbits 
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A complete formalism for constructing initial data representing black-hole binaries in quasi- 
equilibrium is developed. Radiation reaction prohibits, in general, true equilibrium binary con- 
figurations. However, when the timescale for orbital decay is much longer than the orbital period, 
a binary can be considered to be in quasi-equilibrium. If each black hole is assumed to be in 
quasi-equilibrium, then a complete set of boundary conditions for all initial data variables can be 
developed. These boundary conditions are applied on the apparent horizon of each black hole, and 
in fact force a specified surface to be an apparent horizon. A global assumption of quasi-equilibrium 
is also used to fix some of the freely specifiable pieces of the initial data and to uniquely fix the 
asymptotic boundary conditions. This formalism should allow for the construction of completely 
general quasi-equilibrium black hole binary initial data. 

PACS numbers: 04.20.-q, 04.25.Dm, 04.70.Bw, 97.80.-d 



I. INTRODUCTION 

The problem of calculating the gravitational wave- 
forms produced during the final plunge and coalescence 
of a pair of black holes currently is receiving consider- 
able attention. This attention is motivated by the im- 
minent deployment of the Laser Interferometric Gravi- 
tational Wave Observatory (LIGO) and other, similar, 
detectors. It is further motivated by the fact that black- 
hole binaries are expected to be strong sources of gravi- 
tational waves and are considered one of the most likely 
candidates for early detection. The detection and inter- 
pretation of black hole mergers will be greatly facilitated 
by theoretical predictions for the gravitational waveforms 
produced by these events. However, in the absence of pre- 
dictions for the full waveforms, even information about 
the behavior of the orbits will be helpful. In particular, 
information about the properties of the so-called last sta- 
ble orbit, or innermost stable circular orbit (ISCO), is of 
interest. 

Currently, large scale numerical simulations are the 
only known approach for computing the gravitational 
waveforms produced during the final plunge and coales- 
cence of black-hole binaries. Any dynamical simulation 
must start with some initial data, and the results of such 
a simulation are entirely dependent on how well the given 
initial data correspond to the desired physical situation. 
Thus, astrophysically realistic initial data are required 
if we are going to simulate the gravitational waveforms 
of astrophysical black-hole binaries. Because radiation 
damping tends to circularize binary orbits, orbits that 
are nearly circular are of greatest interest. Of course, en- 
ergy loss by gravitational radiation also means that no 
truly circular orbits can exist for black-hole binary sys- 
tems. To keep this fact in mind, we will always refer to 
such orbits as being quasi-circular. 
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The usefulness of current numerical simulations of gen- 
eral black-hole binary systems is limited by two impor- 
tant facts. First, the simulations themselves are only 
beginning to reach the stage where they can produce 
useful information. Furthermore, the accuracy of this 
information is still questionable. Second, the initial data 
currently being used for these simulations is rather crude. 
Initial data for binary black hole systems have been avail- 
able for some time now [fll ||. The approaches for gen- 
erating these initial data can produce rather arbitrary 
configurations and a way of picking which data sets repre- 
sent binaries in quasi-circular orbits was needed. An ap- 
proximation scheme for locating quasi-circular orbits was 
developed H . Based on an effective-potential method, 
this scheme was first applied to non-spinning, equal-mass 
binaries using the inversion-symmetric data of Ref. [Q. 
Later, this approach was applied to non-spinning, equal- 
mass binaries based on puncture dataQ and to spinning, 
equal-mass black holes Q. 

Unfortunately, there are problems with these data sets. 
Because the schemes used to generate this data assume 
the metric to be conformally flat, they cannot produce 
astrophysically accurate data for spinning holes Q. Fur- 
thermore, the estimates for the location of the ISCO 
produced by these schemes do not agree well with re- 
sults based on various post-Newtonian approximation 
methods H 0. There is some speculation that this dis- 
crepancy is rooted in the approximation of conformal flat- 
ness. While the conformal flatness approximation cer- 
tainly introduces some error, I suspect that the most 
important source of error comes from a failure to accu- 
rately approximate the masses of the black holes as they 
approach the ISCO. 

Initial-data schemes relying on the conformal flatness 
approximation have been popular because they are rela- 
tively simple and inexpensive from a computational point 
of view. When maximal slicing is used, an analytic so- 
lution of the momentum constraints exists^, E|. This 
means that only the Hamiltonian constraint, a single 
quasi-linear elliptic PDE, need be solved. Recently, ef- 
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forts have been made to move beyond the approximation 
of conformal flatness (cf. Rcf. (To[ 0). Unfortunately, 
these data have not yet been used to estimate the loca- 
tion of the ISCO. 

For the case of neutron-star binaries, a very different 
approach has been used to construct initial data rep- 
resenting neutron stars in quasi-circular orbits and to 
locate the ISCO. At the core of this approach is the 
assumption that the space-time admits an approximate 
Killing vector, I = (d/dt) + £l(d/d(f>). This approach 
was first tested for the case of a single spinning neu- 
tron star[^2|. Subsequently, variations of this method 
were used to study neutron-star binaries in quasi-circular 

orbits H 0, [l| [tf| © EH- This a PP roach for 
generating initial data for neutron-star binaries in quasi- 
circular orbits is particularly attractive when compared 
to the effective potential method used for black-hole bi- 
naries. First, if the neutron-star matter is to be in hydro- 
static equilibrium, then the 4- velocity of the matter must 
be proportional to the approximate Killing vector, £. The 
assumption that £ is a Killing vector ultimately yields a 
relativistic Bernoulli equation that governs the behav- 
ior of the matter and, most importantly, yields a unique 
value of the angular velocity, f2, of the binary system in a 
circular orbit at a given separation. The presence of mat- 
ter also sets a natural length scale for the the problem 
and conservation of baryon number makes it straightfor- 
ward to connect initial-data solutions at different radii to 
produce meaningful evolutionary sequences. 

(hereafter GGB1 



Recently, Gourgoulhon et al.\2\ 



and GGB2 respectively) have attempted to adapt the 
approach used to construct neutron-star binary initial 
data to the case of black-hole binaries. Their initial re- 
sults are very encouraging, showing much better agree- 
ment with post-Newtonian techniques for locating the 
ISCO. However, there seem to be problems with their 
solutions. It seems likely that their initial data represent 
a much better approximation to astrophysically realistic 
black-hole binaries in quasi-circular orbit than the data 
produced using the effective potential technique. How- 
ever, these new solutions do not represent valid general 
relativistic initial data because they do not satisfy the 
constraint equations everywhere. Furthermore, as I will 
justify later, I do not believe that the particular problem 
being solved in GGB1,2 admits a solution. Interestingly, 
this author had attempted to solve virtually the identical 
problem j2^] but was unable to obtain convergence for the 
coupled system of equations. As the authors point out in 
GGB2, they must regularize their solution for the shift 
in order to obtain an extrinsic curvature that does not 
diverge at the throats of the black holes. It is this reg- 
ularization technique that prohibits their solutions from 
satisfying the constraints everywhere. Although I have 
concerns regarding the validity of the particular solutions 
found in GGB2, I believe that the authors of GGB1 may 
have made a significant contribution to the effort to con- 
struct astrophysically realistic black-hole binary initial 
data. This contribution is rooted in their approach for 



determining the orbital angular velocity, fi. 



The problem with the approach taken in GGB1 is that 
it is not directly applicable outside the restrictions of con- 
formal flatness and maximal slicing, and it seems that a 
solution within these limits does not exist. The limita- 
tions in the GGB1 method are rooted in the boundary 
conditions posed on the inner boundaries of their solu- 
tion domain, that is, the black-hole boundaries. In GGB1 
the authors choose to impose a particular isometry condi- 
tion on the spatial hypersurface. This isometry condition 
naturally yields boundary conditions on the throats of 
the black holes, however, this isometry condition and the 
boundary conditions derived from it have no direct con- 
nection to the desired quasi-equilibrium state of the so- 
lution. Interestingly, the authors of GGB1 choose to mix 
the isometry-derived boundary conditions with a Killing- 
horizon condition that is essentially an equilibrium con- 
dition on the black hole horizon. However, this mixed set 
of boundary conditions is not consistent with the isom- 
etry condition unless the artificial regularization of the 
shift is performed. Again, this regularization prohibits 
the solutions from satisfying the constraints everywhere. 



In this paper, I will develop a formalism for construct- 
ing completely general quasi-equilibrium black hole bi- 
nary initial data. This formalism is also based on a gener- 
alization of the approach used for constructing neutron- 
star binary initial data. The key to this formalism is 
a set of boundary conditions that are based on the ap- 
proximation that the black holes are in quasi-equilibrium. 
Boundary conditions for all of the quantities required to 
construct quasi-equilibrium data can be found based on 
the quasi-equilibrium assumption. These boundary con- 
ditions are completely general. In particular, they do not 
depend on the often-used assumptions of conformal flat- 
ness and maximal slicing. This formalism is completed 
by using the approach developed by GGB1 to determine 
the orbital angular velocity fl of the system, an approach 
that is also rooted in the assumption of quasi-equilibrium. 



In the following sections, I will derive and describe this 
new formalism. Where appropriate, I include compar- 
isons with related methods in order to clarify the ad- 
vantages of the new formalism. In § II, I begin with 
a general discussion of the decomposition of Einstein's 
equations and how the initial-data equations are posed. 
In § [II, I discuss two particular decompositions of the 
constraint equations. In 



[Vj, I discuss the strengths and 



weaknesses of the approach used in GGB1,2 and derive 
the boundary conditions required for the new formalism. 
In § M, I work through the remaining issues involved in 
constructing quasi-equilibrium initial data. And in § VI , 



I conclude with discussion that focuses on the impact of 
possible choices for the conformal 3-geometry and slicing 
condition. 
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II. THE INITIAL-DATA EQUATIONS 

In this work, we will use the standard 3 + 1 decompo- 
sition with the interval written as 

ds 2 = -a 2 dt 2 + j ij (dx i + tfdtydx 1 + (Pdt). (1) 

where 7^ is the 3-metric induced ona( = const, spatial 
hypersurface, a is the lapse function, and (3 l the shift 
vector. The extrinsic curvature of the spatial slice, Kij, 
is defined by 1 



Ki 



(2) 



where C n denotes the Lie derivative along the unit 



normal to the spatial slice, 



Einstein's equations, 



G^v = SitT^l,, then reduce to four sets of equations. 
Two are evolution equations for the spatial metric and 
extrinsic curvature: 



-2aKu + Vi/?j + V,-# 



and 
d t K, 



a[R 



2K i ,K t A 



SirSij + 47T7ij(S - p)] 



(3) 



(4) 



■KuVjf + KjtViP 1 . 
The remaining two are the constraint equations: 

R + K 2 - A , ,/\' ' = 167rp (5) 

and 



(6) 



of the 12 degrees of freedom found in the initial data, 7^ 
and K~ij. The covariant nature of Einstein's equations 
allows us four gauge degrees of freedom. The remain- 
ing four degrees of freedom are the dynamical degrees of 
freedom and are split two apiece between the metric and 
extrinsic curvature. 

The task of constructing initial data for a Cauchy evo- 
lution via Einstein's equations can be broken down into 
a few fundamental task shared by any approach. 

1. A decomposition of the constraints must be chosen 
that specifies how the constrained, gauge, and dy- 
namical degrees of freedom are associated with the 
initial data. 

2. Boundary conditions must be chosen for the con- 
strained degrees of freedom such that a) the con- 
straint equations reduce to a well-posed set of el- 
liptic PDEs and b) they are compatible with the 
desired physical content of the initial data. 

3. A choice for the spatial and temporal gauge must be 
made that is compatible with the desired physical 
content of the initial data. 

4. The dynamical degrees of freedom must be specified 
so as to yield the desired physical content in the 
initial data. 

Once these steps have been completed, all that remains 
are the technical aspects of solving the coupled set of 
elliptic PDEs. The physics lies entirely in the four enu- 
merated steps given above. 



Here, Vj, Rij, and R are, respectively, the covariant 
derivative, Ricci tensor, and Ricci scalar associated with 
the spatial metric 7^ . K = K\, $y is the matter stress 
tensor, S = SI, p is the matter energy density, and j l is 
the matter momentum density. These matter terms are 
related to the stress-energy tensor via: 



3 1* 
P 

T 

-L fl U 



= l^T a0 , (7) 

= - 7 > Q r ra , (8) 

= n"n%„ (9) 

= + 2n( A J„) + n^n v p. (10) 



Initial data for a Cauchy evolution via Einstein's equa- 
tions consists of specifying 7^ and Kij on an initial hy- 
persurface. The Hamiltonian constraint (||) and the mo- 
mentum constraints (^J) constitute the initial-data equa- 
tions of general relativity, and the initial data must sat- 
isfy these equations. Equations (|^) and (^) constrain four 



III. DECOMPOSITIONS OF THE 
CONSTRAINTS 

In the weak field limit where Einstein's equations can 
be linearized, there are clear ways to determine which 
components of the initial data are dynamic, which are 
constrained, and which are gauge. However, in the full 
nonlinear theory, there is no unique way to perform 
the decomposition. Perhaps the most widely used class 
of constraint decompositions are the York-Lichncrowicz 
conformal decompositions which are based on a con- 
formal decomposition of the metric and various other 
quantities pi |25l pel. 



The metric is decomposed into a conformal factor, ip, 
multiplying a conformal metric: 

lij=4>%j- (11) 
Using (pr|), the Hamiltonian constraint (||) can be written 



1 In this paper, latin indices denote .^-dim ension al spatial indices. 
Greek indices are 4-dimensional. In § IV C, latin indices can 



be replaced with indices varying over the two dimensions of the 
surface whenever only 2-dimensional tensors are involved in the 
expression. 



VV - \^R - \^K 2 + ^ 5 K i:j K ij = -2^/> 5 p, (12) 

where V 2 = V J Vi, and Vj and R are the covariant 
derivative and Ricci scalar associated with 7^. With 
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this conformal decomposition, the conformal metric, -yy 
is considered to encode the dynamical degrees of free- 
dom of the metric and the initial gauge choices for the 
three spatial coordinates. The conformal factor, tp 1 is the 
constrained portion of the metric. 

Most decompositions of the extrinsic curvature begin 
by splitting it into its trace, K and tracefree parts, Aij, 



Ki 



A,. 



(13) 



The decomposition proceeds by using the fact that we 
can covariantly split any symmetric tracefree tensor, S lJ , 
as follows [B7|: 



S 13 = (LX) i3 + T l - 



(14) 



Here, T %3 is a symmetric, transverse-traceless tensor (i.e., 
VjT 1 - 7 = and = 0) and, in three dimensions, 



(JLX) tj = S/ i X j + V 3 X l - h lj V e X £ 



(15) 



At this point, the decomposition can proceed in several 
different directions, and we will consider two of them be- 
low. However, in all cases, K is considered to encode 
the initial temporal gauge choice, T %3 encodes the two 
dynamical degrees of freedom, and X 1 are the three con- 
strained degrees of freedom. 



A. Conformal transverse-traceless decomposition 

Historically, the most widely used decomposition ex- 
presses the extrinsic curvature as 



K ij = f 



-10 



(16) 



In this case, the transverse-traceless decomposition takes 
place with respect to the conformal metric. It will also 
be convenient to define a conformal trace-free extrinsic 
curvature, A 13 , by 



A ij ee ^ 10 A 13 = (Lxf 3 + Q lj . 
Using ( |l6| ) in the momentum constraint (^J) yields 



where 



A h X l = l^ l K + 8ml; 10 f, 



A L X* ee v^lxy 3 ' 

= W 2 X l + |V*(VjJC*) + R)x j , 



(17) 



(18) 



(19) 



and Bh is the Ricci tensor associated with 7^ . In terms 
of (|17|), the Hamiltonian constrain (|l2|) becomes 



+ \i,- 7 A l3 A^ 



(20) 



Together, Eqns. QIOJ) and (|16|) (g0|) constitute the full 
conformal transverse-traceless decomposition of the con- 
straints. To solve this set of equations, one must specify 
Jij j Q y j an( i K m addition to the matter terms p and j l . 

The reason that this decomposition has been so widely 
used stems from the simplifications that occur when the 
freely specifiable data are chosen so that the initial data 
hypersurface is maximal, K = 0, conformally flat so that 
jij is a flat metric, and Q lJ is chosen to vanish. Then, 
for vacuum spacetimes, the constraints reduce to 



and 



V 2 ip+±i>- r A ij A i i =0, 



V 2 JT + 3V*(V 3 -Jr J ') = 0, 



(21) 



(22) 



and Vj is, in this case, the usual flat space covariant 
derivative operator. What is remarkable about this sim- 
plified system of equations is that Eq. (^2|) decouples 
completely from (|2l]) and has analytic solutions that rep- 
resent a single black hole with any desired value of linear 
and angular momentum|s , 28 . Furthermore, since (|22| ) is 
linear, any number of these solutions may be superposed 
to represent multiple black holes. All that is required 
to construct complete initial data is to solve the single 
quasi-linear elliptic PDE in Eq. (|2l|). 

All of the work to date that has used the effective po- 
tential method to locate quasi-circular orbits in black- 
hole binary systems^, ||, |) has used initial data com- 
puted by this system. The only difference between the 
data used in Refs. [[| 1|, an d in Rcf. || is in the topology 
of the initial data hypersurfaces. 



B. Conformal thin-sandwich decomposition 

An alternative decomposition recently proposed by 
York [B9| is a generalization of an approach used first by 
WilsoriQ and later by others @, |l|, |l| for the study 
of neutron-star binaries. Furthermore, aside from some 
trivial differences in conformal scalings, it also subsumes 
the approach outlined in GGB1. This decomposition dif- 
fers from the conformal transverse-traceless decomposi- 
tion significantly in philosophy, but is very similar in 
form. The difference is found in the decomposition of 
the extrinsic curvature. In the conformal thin-sandwich 
decomposition, Ku is decomposed as 



K t3 = 



or 



, V| ((UT-u^+^K, (23) 



A 13 = 



(24) 



Here, u %3 takes the place of Q 13 in representing the freely 
specifiable dynamical portion of the extrinsic curvature. 
However, owing to the nature of the thin sandwich de- 
composition, it also has a more physical meaning as we 
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shall see in a moment. Notice that in Eq. fl24), it is not 
the conformal trace-free extrinsic curvature that is de- 
composed, but rather 2aA l K The new scalar, a, is the 
conformal lapse defined by 



(25) 



Also notice that, in the transverse-traceless decomposi- 
tion, the vector X 1 has been replaced by the shift vector, 

Using the thin-sandwich decomposition ( p^ ) in the mo- 
mentum constraint flra) yields 



±aip 6 V l K (26) 
+ aVj(lu ij ) +16 7 r5V 1 V- 



Taken together, Eqns. (|ll|), (|23|), (]2J), (|2J), (|19|), and 
( pp| ) constitute the full thin-sandwich decomposition of 
the constraints. To solve this set of equations, one must 
specify 7^ , u lJ , K , awe? a in addition to the matter terms 
p and j\ 

The thin-sandwich decomposition differs from most 
initial-data decompositions in that its derivation uses a 
second spatial hypcrsurface that is infinitesimally sepa- 
rated from the primary spatial hypersurface. The time 
and spatial coordinates in these two hypersurfaces are 
connected by the kinematical variables, a and j3 l , and 
this connection allows us to obtain a physical interpre- 
tation for u lJ , the freely specifiable dynamical portion of 



the extrinsic curvature. The result is that 



is a sym- 



metric, tracefree tensor with the physical interpretation 
that 



dtjij, 



(27) 



and t is the time coordinate associated with the time-like 
vector f that connects points in the two adjacent hyper- 
surfaces that carry the same spatial coordinate labels, 
i.e. 



(28) 



and n M is the timelike unit normal to the primary spatial 
hypersurface. 

This decomposition has proven useful for construct- 
ing quasi-equilibrium initial data sets because, in the 
frame where the configuration is assumed to be time- 
independent (the corotating frame), it is a natural ap- 
proximation to set Uij = 0. To my knowledge, all of 
the work so far (including Refs. f| ||, |T|, |o[ |f| , and 
GGB1,2) use the same set of assumptions to fix the freely 
specifiable data in the thin-sandwich decomposition. In 
particular, they all choose the initial data hypersurface 
to be maximal, K = 0, conformally flat so that 7^ is a 
flat metric, «, 3 is chosen to be zero, and they all fix the 
lapse by demanding that dtK = 0. The final assumption 
is a reasonable one since a quasi-equilibrium solution was 
the goal of each of these works and the time vector ( |2"§| ) 
is the one corresponding to the corotating frame. 



This conformal thin-sandwich decomposition will form 
the foundation of the initial-data formalism being devel- 
oped. However, unlike previous works, we will only make 
the restriction that Cty = 0. In particular, we will not re- 
strict the formalism to conformally flat maximal slices of 
the full spacetime, and we will require that all boundary 
conditions and any further restrictions apply (at least in 
principle) to general slicings and conformal 3-geometries. 



IV. BOUNDARY CONDITIONS FOR THE 
CONSTRAINT EQUATIONS 

A. The GGB approach 

We now return to the specific problem of using the 
thin-sandwich decomposition to construct initial data for 
a pair of black holes in quasi-circular orbits. In addi- 
tion to the standard assumptions outlined at the end 
of § IIIB, GGB1 also makes the assumption that the 

That 



initial-data hypersurface is inversion symmetric, 
is, the initial data hypersurface is taken to consist of two 
asymptotically flat hypersurfaces that are joined together 
at the throats of each black hole, and are isometric to 
each other with each throat being a fixed point set of the 
isometry condition. The isometry condition then imposes 
a set of boundary conditions at the throats on all quan- 
tities on the manifold. There are actually two choices 
for the boundary conditions, and GGBI make the choice 
that requires the lapse to vanish on the throats. For the 
shift vector, the isometry imposes the following boundary 
conditions 



d(3 e 

d r 



0, 
0, 

0. 



(29) 
(30) 

(31) 



where we are assuming the throat is a coordinate sphere 
and |s denotes evaluation on the throat. 

Since this version of inversion symmetry requires that 
a I s = 0, regularity of the extrinsic curvature at the 
throats requires 



(L/3H 5 = 0. 



(32) 



As GGBI point out, this condition is potentially prob- 
lematical. 

Given the choice for the isometry condition used in 
GGBI, it follows |3^] that the throats are apparent hori- 
zons. GGBI make the interesting choice of using this 
fact to change boundary conditions (|3^) and ([5l]) on the 
shift. They obtain new boundary conditions on the shift 
by demanding that d/dt be null on the throats so that 
the throats are Killing horizons. For the case that the 
lapse vanishes on the throats, this condition yields 



/9*|s = 0. 



(33) 
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Altering the boundary conditions on the shift in this way 
poses a problem. Using Eq. ( {j3| ) to replace conditions 
( |30| ) and means that the solution is not guaranteed 
to be inversion-symmetric. This, in turn, means that the 
inner boundaries are not guaranteed to be apparent hori- 
zons, from which it follows the Killing horizon boundary 
conditions may not be appropriate for the inner bound- 
aries. 

In GGB1, the authors go to great length to show that 
inversion-symmetry, their Killing-horizon boundary con- 
ditions on the shift, and the required regularity condition 
( |32| ) can be compatible. However, satisfying the Killing- 
horizon boundary conditions together with some of the 
inversion-symmetry boundary conditions does not guar- 
antee that the regularity condition ( |32| ) is satisfied, or 
that the resulting solution is inversion symmetric. If fact, 
as pointed out in GGB2, it seems that their solutions are 
not compatible with this condition without some kind 
of artificial regularization of their solution for the shift. 
This is not surprising. We can rewrite the evolution equa- 
tion for the metric as 



(34) 



Uij + §7ij (V fc /3 fe + 6f3 k V k lm[>- ^ & aK 



j,y = tp^Uij. Given the the Killing- horizon bound- 



with Ui 

ary conditions and the choices of maximal slicing and 



0, Eq. fl3J) reduces to 



d t jij\ 



2 ,4- d ^ 



(35) 



It is demonstrated in GGB1 that regularity requires 
d(3 r /d r \s = and the authors take Eq. (35) as estab- 



lishing regularity. However, the choice Uj 



is not 



equivalent to dtfij = 0, but only fixes the time rate of 
change of the conformal metric. The time rate of change 
of the conformal factor is not fixed by this condition and 
is only guaranteed to vanish if t is a true Killing vector. 
Thus, for a quasi-equilibrium solution of the constraints, 
it seems reasonable that df3 r /d r \s is small in some ap- 
propriate norm, but it is unlikely to vanish. 

So, the approach outlined in GGB1 has two techni- 
cal problems: the boundary conditions being used do 
not guarantee that the solution is inversion symmetric or 
that the solution will be regular at the inner boundary. 
Both of these problems are fixed in GGB2 by an artifi- 
cial regularization of the shift solution. The unfortunate 
side effect of this regularization procedure is that the 
solutions no longer satisfy the constraint equations and, 
thus, cannot represent valid initial data. In fact, it seems 
that a regular, inversion-symmetric solution for the ex- 
trinsic curvature is not possible (or is at least extremely 
unlikely), given the assumptions of GGB1. 



B. Killing- horizon boundary conditions 

Although the particular set of boundary conditions 
proposed in GGB1 does not seem to yield a well-posed 
system of elliptic equations, the idea of using some kind 
of boundary conditions on the apparent horizons of the 
black holes is promising. An alternate approach for ob- 
taining boundary conditions on the shift is found by a 
slight generalization of the the Killing-horizon conditions 
proposed in GGB1. As we will see, this change in per- 
spective will allow us to derive boundary conditions on 
the shift that will yield corotating black holes, irrota- 
tional black holes, or black holes with some arbitrary 
angular velocity. We begin by noting that the outward 
pointing null vector whose vanishing expansion defines 
the apparent horizon can be written as 



(36) 



Here, s M is the outward pointing unit vector normal to 



the apparent horizon (see § IV C for further details). In 
terms of components, this can be written as 



ft" = 



(37) 



where I have chosen a normalization so that fc* = 1 . Note 
that is null by construction and k^k^ — for any 
choice of the shift vector. 

If we assume rotation in the 4> direction, then the ap- 
proximate killing vector for a quasi-equilibrium configu- 
ration is given by 



(38) 



where f2 is the orbital angular velocity of the binary, and 
d/dcj) represent the asymptotic Killing vector for rotation 
in the <f> direction. In the corotating coordinate system, i 
is the time vector and it has components £^ = [1, 0, 0, 0]. 
If the individual black holes are corotating with the sys- 
tem, then they do not appear to be rotating in that coor- 
dinate system. In this case, the Killing-horizon condition 
implies that £^ should be parallel to /c M which immedi- 
ately yields 



as |s. 



(39) 



Notice that this condition reduces directly to Eq. (|32 
when a\s = 0. 

We can also derive a boundary condition for black holes 
that are not rotating with respect to observers at rest 
at infinity, black holes that we might loosely refer to as 
being irrotational. In this case, the outward pointing 
null vectors at the apparent horizon should be parallel 
to t, not to £. In the corotating coordinate system, the 
components of t can be written symbolically as 



f = [i,-n(d/d(f>y 



(40) 
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It follows immediately that the appropriate boundary 
conditions for irrotational black holes are 



(41) 



For both boundary conditions, corotation given by 
(|39|), and irrotational given by (j4l|), the conditions are 
given in terms of the corotating coordinate system. This 
means that the appropriate asymptotic boundary condi- 
tions for the shift arc that 



0%- 



(42) 



The corotation and irrotational conditions can be gen- 
eralized to yield a boundary condition for an equilibrium 
black holes with arbitrary spin. Let \ l denote a flat-space 
Killing vector for rotation in an arbitrary direction, just 
as (d/d(f>) 1 represents rotation about in the <j> direction. 
With Q, x representing the angular velocity of rotation of 
the black holes as measured by observers at rest at in- 
finity, we find that the components of the null vectors 
generating the horizon are 



(43) 



when written in terms of the corotating coordinate sys- 
tem. This immediately yields the boundary conditions 
for a quasi-equilibrium black hole with arbitrary spin: 



P\s=aJ\ s -n x x t \s + n [4r 



(44) 



As we will see below, we will need to be a bit more care- 
ful in defining the shift boundary condition for the irro- 
tational and general spin cases, but (|4l]) and ( fii]) give us 
a good starting point for these conditions. 



C. Quasi-equilibrium boundary conditions 

The choice of the particular inversion-symmetric class 
of black- hole initial data used in GGB1 was made because 
it guarantees that the inner boundary surfaces will be ap- 
parent horizons. If we are to give up inversion-symmetry, 
we must find some other way of fixing boundary condi- 
tions at the inner-boundary surfaces for all of the quan- 
tities needed to construct initial data. These boundary 
conditions can be found by forcing the inner boundary to 
be an apparent horizon and by assuming that the associ- 
ated black hole is in quasi-equilibrium. Pieces of the re- 
sulting approach have been worked out by Thornburg [[33| 
and by EardleyQ, but to my knowledge have not been 
adapted for constructing quasi-equilibrium initial data. 
Because there are several errors in the derivations and 
results in Ref. p4| , and to provide a uniform notation, I 
will rederive all of the equations here. 



Our inner boundary surface, S, is assumed to be a 
spacelike 2-surface with topology S 2 . Because S is closed, 
we can define its unit normal via 



Si = AV,t, 



(45) 



where A is a normalization constant fixed by s l Si = 1 
and t is a scalar function which has a level surface t = 
that defines the surface S. The 4-dimensional general- 
ization of s l has components s M = [0, s l ] obtained from 
the condition that s^n^ = 0. 

The metric, hy, induced on S by 7y is given by 

hj = Jij - SiSj. (46) 

We also define the extrinsic curvature, Hij , of S embed- 
ded in the 3-dimensional spatial hypersurface as 

Hij = -h>;hp {kSe) = -\C.htj. (47) 

Naturally associated with S are two sets of null vectors: 
a set of outgoing null rays, and a set if ingoing null 
rays, fc M , defined by 

fc' 1 = -2= (r^ + s^) and U l = (n* - s^) . (48) 



V2 



Associated with each set of null rays is an extrinsic curva- 
ture of S as embedded in the full 4-dimensional manifold. 
These are defined as 



(49) 
(50) 



To simplify the following, we will introduce various pro- 
jections of Kij along and normal to S: 



Jij = h^h l jK k i, 



kit 



(51) 
(52) 
(53) 

= J i:j + 2s {i Jj ) + SiSj(K - J). (54) 
We can then simplify Eqns. (^) and (^) to 

Now, wc define the expansion of outgoing null rays, <r, 
and ingoing null rays, &, via 



a = h^ l3 = ^(J + H). 
6 = h ij ±ij = -^{J-H). 



(56) 
(57) 



Finally, we will define the shear of the outgoing null rays, 
Oij , and ingoing null rays, &ij, via 



hha<j and a 



\hij&. (58) 



In order to generate boundary conditions, we need to 
make some assumptions. First, we will demand that our 
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inner boundary, S, is an apparent horizon. This is equiv- 
alent to demanding that the expansion of the outgoing 
null rays vanishes on S. Thus our first condition is that 

a = 0. (59) 

If we are looking for quasi-equilibrium configurations, 
then we want S to remain at the same coordinate loca- 
tion in the 3-dimensional hypersurface as time passes. To 
enforce this condition, it is necessary that Cfj = dtT = 0, 
where t is the approximate Killing vector associated with 
our demand of quasi-equilibrium. However, there is no 
reason that coordinates cannot be free to move on S. In 
fact, this freedom is necessary to allow for rotation of the 
black hole. This means that, in order to keep S at the 
same coordinate location, dtT = is too strong a condi- 
tion. If we define £ as the part of t that is orthogonal to 
5, then it is necessary and sufficient that 



£ c t = and HViCcr = ViCtT = 



(60) 



where Vj is the 2-dimensional covariant derivative com- 
patible with hij. If we define the normal component of 
the shift as 

#L = (61) 

then we can write ( as 

C 1 = q»/ + /W, (62) 

We now need to consider how the expansions, a and 
cj, evolve along £. For this calculation, a few identities 
are crucial: 

h ij Rij = !(.R- H 2 + H i3 H i:i + R), (63) 
s V'i?y = \{R + H 2 - H i3 H i3 -R), (64) 
A[£ c ,V M ]r = -fi.KVA+fctJf-/)], (65) 

where R is the Ricci scalar associated with and 
is the 4-dimensional covariant derivative. The first two 
identities are obtained from various combinations and 
contractions of the Gauss-Codazzi-Ricci equations that 
govern the embedding of S in the the 3-dimensional hy- 
persurface. The third identity can be obtained by direct 
calculation. 

As useful intermediate results, we obtain: 



-Jif + 8irS ij s i s j + |R - V' 
J lJ H i3 +JH- HK - 8tt j t s l 



-Vi J 1 + 2J J V, 



p ± +Hs i V i a, (66) 



C C H = 



J tj H %i - 8nj iS z - VjJ 8 - 2J l Vi a 
|(J 2 + H 2 ) + ^(JijJ i3 + H t jH i3 



J t J l -JK + 87T/5 - |R + V 2 



+Js l ViOi. 



(67) 



These equations follow from the identities (|6^)-(|65|) and 
the evolution and constraint equations Finally, 
we can simplify these to the desired results: 



£c<7 



l 

1 

73 



1 
1 

V2 



°{° + \6- -feK) + £\ [fi ± + a) 



+8wT^k' 1 k v {fix - a) + (Ts'Via, (68) 
6(6 + \<r- zj^K) + £~\ (/? x - a) 
6(¥-\°--\K)+V 

+8irT fiU k fl k ,/ ] (Px + a)~ as^ia, (69) 



where 



V 
V 
£ 



W(V l + J l )(V 3 +J J ) 

hVfyi-JiXVj-jj) 



2 R i 

l-R, 



(70) 
(71) 
(72) 
(73) 



We note that £ and £ are both non-negative so long as 
the matter satisfies either the weak or the strong energy 
condition. The two equations for the evolution of the 
expansions, (|6^) and (|6^), are not yet restricted by any 
assumption of quasi-equilibrium. The only assumption in 
their derivation is that the surface S remain at a constant 
coordinate location which is embodied in ( |gc| ) . 

If we now restrict Eqns. (68) and ( |69| ) to apply to ap- 
parent horizons by enforcing condition (^9|) and also im- 
pose quasi-equilibrium on the evolution equations for the 
expansions, 



C^o = and Cc6 = 



(74) 



then we find 

= 

&s l Via = 



V(p ± - a), 

6(6--kK) + & ij 6 i - 



(75) 



i 

71 



(fi± - «) 



1 

V2 



a(\a-^K)+V (/? 



(76) 



These equations are based on the approximation that 
the shear of the outgoing null rays on the apparent hori- 
zon vanishes, which is true for the case of a stationary 
black hole, and on the requirements that T^k^W — 
T^k^ky — TavWh" — on the apparent horizon if the 
black hole is in quasi-equilibrium. Finally, I reempha- 
size that the conditions for quasi-equilibrium on the two 
expansions is given by Eq. ([74]), not by dt<J = dt& = 0, 
so that the location of the horizon is fixed but the black 
hole is allowed to rotate in S. 

On inspection, Eq. (f75|) yields the solution 



a, 



(77) 
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which is compatible with the Killing horizon boundary 
condition obtained in the case of corotation (|39|). The 
Killing horizon boundary conditions for irrotation and 
general rotation, (41) and (Q), are not compatible with 
([Ft]) unless we restrict the rotational terms to act in the 
surface S. Thus, the correct boundary condition on the 
shift for an irrotational black hole should be 



P i \s = as i \s + Qhi ( — 



8 v 

00 



and the condition for general rotation should be 



as Is 



o,x 3 - n 



d 



(78) 



(79) 



Notice that our assumption of quasi-equilibrium has lead 
to one condition on the shift, and yet Eqns. (|3^), (|7Sl) 
and ( |79| ) specify three conditions on the shift. This is 
appropriate since we have made specific choices for the 
rotation of the black hole, namely corotation, irrotation, 
or some general rotation specified by tt x and \ l - 

The condition that S be an apparent horizon, (59), 



also yields a useful boundary condition if we work in the 
conformal space. The conformal transformation on Yy 
( pi] ) induces a natural conformal weighting for hij and 
for the unit normal to S, 



'ij 
s i 



ip hij, 



We find that Eq. (p9[) reduces to 



1 j 



■1?J)\t 



(80) 
(81) 

(82) 



which is a non-linear boundary condition that incorpo- 
rates both ip and its normal derivative. While this bound- 
ary conditions looks somewhat complicated, it reduces to 
the often used minimal surface boundary condition when- 
ever J\s = 0, and it has previously been used successfully 
to construction initial-data sets [p3| . 

We can now further simplify the remaining condition 
given in Eq. (f76[). From Eqns. (p 7 ^), (|57|), and (|59|), we 
obtain 

a = V2J. (83) 
(§(]), and @, this reduces @ to 
= -V 2 (J 2 - JK + V)a\ s , (84) 



QUASI-CIRCULAR ORBITS FOR 
BLACK-HOLE BINARIES 



The primary goal of this paper is to fully define a for- 
malism for constructing initial data sets representing as- 
trophysically realistic black- hole b inaries in quasi-circular 
orbits. As was shown in 8 IIIB, the thin-sandwich de- 



composition of the constraint equations appears to be 
a natural choice to use for this purpose. In particular, 
the simple choice of u 13 = is r equire d by the assump- 
tion of quasi-equilibrium. In § [V C , boundary condi- 
tions compatible with quasi-equilibrium were developed. 
However, there are still several freely specifiable quanti- 
ties that we have not considered within the context of 
quasi-equilibrium. In particular, we must consider how 
the initial temporal gauge choice evolves off of the initial 
slice via a. We must decide what value to use for the or- 
bital angular velocity, in setting the outer-boundary 
conditions ([l2|). And finally, we must choose the confor- 
mal 3-geometry, 7y, and the initial temporal gauge via 
K. 



A. The temporal gauge choice 

In all previous work on quasi-equilibrium data for ei- 
ther neutron stars or black holes, maximal slicing has 
been used. With this choice, K — 0, and it is natural to 
use the "maximal slicing equation" to fix the lapse. The 
maximal slicing equation has been used extensively in 
numerical evolutions and is simply a linear second-order 
elliptic equation for the lapse derived from the evolution 
equation for K by setting dtK — (cf. Ref. p5|). For a 
quasi-equilibrium situation, it is natural to use the condi- 
tion dtK = 0, but to generalize it to any arbitrary value 
for K. In this case, the equation that results, 



V 2 (cn/0 - a 



\i>R + ^ip 5 K 2 + l^- 7 A l3 A 13 (86) 
+2irip 5 K(p + 25)] = ^ 5 f3 l ViK, 



Together with (jJJ 
with 



is often called the "constant trace- if equation" (constant 
refers to the fact that K is constant in time, not constant 
in space). 

For a quasi-equilibrium situation, we would like to have 
dtlij — and d t Kij — 0. However, this will only be ob- 
tainable in true equilibrium situations where we have ex- 
act Killing vectors. We can decompose these conditions 
to get dt'jij = 0, dti> = 0, d t Aij — and d t K = as the 
possible quasi-equilibrium conditions that we might ap- 
2V 2 ln-0]. (85) ply- Of these, only d t jij = Uij is part of the freely speci- 
fiable data of the thin-sandwich decomposition. Given 
Here, Vj and R are the covariant derivative and Ricci that we must fix 7^ and that we choose Uij — as a 
scalar associated with the conformal metric hij. Equa- 



tion (|84j) represents a complicated mixed condition on 
the lapse that involves an elliptic operator acting over 
the closed surface iS. There is no guarantee that ( pH ) is 
invertible, however, if it is, then (B4j) represents a viable 
boundary condition. 



quasi-equilibrium condition, it is not reasonable to as- 
sume that we can also find dtAij = except in true 
equilibrium conditions. In addition to two dynamical de- 
grees of freedom, 7^ also encodes the initial spatial gauge 
choice. Similarly, K encodes the initial temporal gauge 
choice and parity with the spatial gauge suggests that we 
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should let quasi-equilibrium fix d t K = and that it is not 
reasonable to expect d t tp = except in true equilibrium 
conditions. 

We can examine the condition dtip = even further. 
From Eq. ([m]) we find that this condition is equivalent 
to 



Vfe/3 fc + 6/3 fe V fc ln^ 



ip 6 aK 



0. 



(87) 



We might consider using ( j87| ) to fix K instead of choos- 
ing it freely. However, if we use (|87j ) to replace K in 
either Eq. (|2^) or (|l8|) we find that the resulting elliptic 
operator is non-invertible. One might then hope that an 
iterative scheme could be used to fix K so that dttp = 
is satisfied. However, numerical experiments with this 
approach have shown it to be unstable |3(|. 

Given that the thin-sandwich decomposition requires 
that a be fixed, it is fortunate that the condition of quasi- 
equilibrium provides an elliptic equation that fixes this 
quantity, just as the condition of quasi-equilibrium pro- 
vides the means of constructing boundary conditions for 
the constraints. Furthermore, quasi-equilibrium also pro- 
vides a boundary condition on a for use in (|S6|). This 
boundary condition is given by Eq. ( |S4| ) with a replaced 
by a via (£51). At spatial infinity, we have that 5 = 1. 



that ip = 1, 5 = 1, and 
is well-posed, solving it wil 



B. Obtaining circular orbits 

Consider solving the thin-sandwich equations, ( pOj ) and 
(26), in conjunction with the constant trace- K equation, 
(36). Assume that we fix a conformal 3-geometry 7^ (not 
necessarily flat), that we fix K (not necessarily a maxi- 
mal slice), and that we choose — 0. We will demand 
that the apparent horizons associated with each hole in 
our black-hole binary system occur on specified closed 2- 
surfaces in the given 3-geometry. On these inner bound- 
aries, we will use the conditions given in Eqns. (|82]), (84), 
and (^) or (|7^) together with the asymptotic conditions 

. Assuming that the system 
result in two one-parameter 
families of solutions, one for the corotation condition ( |39| ) 
and one for the irrotational condition (fT8|) . These fami- 
lies of solutions are parameterized by the orbital angular 
velocity, fi. 

As is pointed out in GGB1, it is unreasonable to as- 
sume that the entire family of solutions will satisfy the 
physical condition of quasi-equilibrium in which we are 
interested. This situation is similar to that encountered 
in searching for quasi-circular orbits via the effective po- 
tential method^]. In fact, the effective potential in that 
approach was chosen on physical ground to pick those so- 
lutions which represented approximately circular orbits. 
However, that effective potential was necessarily some- 
what ad-hoc. In my opinion, the most important con- 
tribution of GGB1 is their more physically well-founded 
and covariant method for determining which members 
of the family of solutions satisfy the condition of quasi- 
equilibrium. 



The condition for choosing appropriate values of CI 
given in GGB1,2 is justified largely in terms of its relation 
to a general relativistic version of the virial theorem fl37|. 
The condition is that Cl is chosen so that the standard 
ADM mass, E ADU , and the Komar(or KVM) massQ, 
M K , agree. In my view, the justification for this condi- 
tion resonates most clearly when one considers that, in 
general, E ADM ^ M K . However, as was shown by Beig 39 
(see also Ref. |Q), these two mass measures agree when 
the spacetime is stationary. Since a quasi-equilibrium 
solution is approximating a stationary spacetime, it is 
natural to demand that E Aml — M K , and in this case, we 
should obtain quasi-circular orbits. 

Since, for a general gauge choice, E ABM is not necessar- 
ily given by the monopole piece of the conformal factor, 
we will use a general definition of the ADM mass Ell, 



e a »u = j^4 -"v,,(t;;' 



^G)d 2 S 3 



(88) 



Here, Qij = 7y — fij, fij is the flat metric to which 7^ 
asymptotes, and d 2 Si is the covariant surface area ele- 
ment. We note that in (|8^), indices can be raised and 
lowered with either 7^ or /y, the trace, Q, can be ob- 
tained with either metric, and Vj can be replaced with 
the flat covariant derivative. Similarly, a general defini- 
tion of the Komar mass can be written as |p7| 

M « = ^ j 7 iJ (V 4 a - (3 k K lk )d 2 S r (89) 

In many cases, f3 k Kik will fall off faster than 0{r~ 2 ) and 
the second term in (89) will not contribute. However, in 
some gauges this term is important. An example is the 
Painleve-GuUstrand coordinate system (cf. Ref. W2[ ) . 



C. Quasi-equilibrium evolutionary sequences 

It is highly desirable to be able to construct evolu- 
tionary sequences of quasi-equilibrium binary configura- 
tions. In particular, this facilitates locating the ISCO. 
When one considers neutron-star configurations, there is 
a natural way to connect neighboring solutions. Namely, 
the number of baryons contained in each star should not 
change as the star secularly evolves to smaller separation. 
With black holes, there is no such conserved quantity and 
one must fix some normalizing condition. 

In order to construct evolutionary sequences from 
quasi-circular orbits obtained by the effective-potential 
method, both the normalizing condition and the orbital 
angular velocity had to be determined by some means. 
One condition that should be applicable to evolutionary 
sequences of quasi-equilibrium configurations |?], [l3| Q is 



Cl 



dE, 



dJ 



(90) 



This relation has been shown to hold well for evolution- 
ary sequences of both corotatingjl4| and irrotational [p0| 
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neutron-star binaries. Because the effective-potential 
method cannot determine the orbital angular velocity 
of each configuration of holes in circular orbit, Eq. j90| ) 
must be used to determine Q. But, locating each circular- 
orbit configuration, as well as determining requires a 
normalizing condition. In work so far with the effective- 
potential method, this normalization condition has been 
chosen to be an ad hoc definition for the mass of an indi- 
vidual black hole in the binary system, with the further 
assumption that this mass remains constant as the black 
holes secularly evolve closer together. While this defini- 
tion had appropriate limiting behavior when the black 
holes have large separation, this normalization condition 
seems problematical in the strong-field limit. 

However, because quasi-equilibrium data constructed 
from the combination of the thin-sandwich and constant 
trace- -R" equations naturally yields a value for £1 for each 
circular-orbit configuration, GGB2 point out that ( |90| ) 
can be used to set the normalization condition for con- 
structing evolutionary sequences. This appears to be 
a much more well-founded condition than the ad hoc 
condition proposed in Ref. || for the effective-potential 
method. That ad hoc condition is essentially rooted in 
the notion that the area of the apparent horizon, an esti- 
mate of the true irreducible mass of a black hole, should 
remain constant during secular evolution. Interestingly, 
the first results found in GGB2 show agreement with the 
ad hoc method to the level of numerical error in the solu- 
tions. Given the proble ms wi th the solutions obtained in 
GGB2, as outlined in § IV A, it is too soon to make any 



for the binary problem does not seem possible when 
inversion-symmetric boundary conditions are chosen be- 
casue a |s = 0. If we consider the case of a single black 
hole, we can easily find a CFMS solution of Eqns. (|2(|), 
(p(f), and ( |86| ) that satisfies the inner boundary condi- 
tions given by Eqns. @, (H), and (§J). This solution 
is simply the Schwarzschild solution in isotropic coordi- 
nates. In this case, although a = on the apparent 
horizon, the solution is regular. This is possible because 
the spacetime admits a true Killing vector. The direct 
generalization of this solution to a binary system yields 
the problem attempted in GGB1,2 that does not yield 
a regular solution. However, this does not mean that a 
quasi-equilibrium solution is impossible given the CFMS 
assumptions. The nonlinearity of the system of equations 
and boundary conditions allow for at least the possibility 
of such solutions, even if it seems unlikely. 

The apparent problems with constructing a quasi- 
equilibrium, binary, CFMS solution should not be taken 
to suggest that conformally flat quasi-equilibrium binary 
solutions are unlikely. The problem with the CFMS as- 
sumptions is with maximal slicing, not conformal flat- 
ness. Again, for the case of a single black hole, we 
can easily find solutions of the quasi-equilibrium equa- 
tions and boundary conditions that are conformally flat. 
One example is obtained from the Schwarzschild solution 
in Kerr-Schild coordinates, also referred to as ingoing 
Eddington-Finkelstein coordinates (cf. Ref. (4^]). With 
the radial coordinate transformation 



conclusions regarding the validity of the ad hoc condition, 
but this point should be considered further. 




^(t-vW), (91) 



VI. DISCUSSION 

So far, we have not considered specific choices for the 
conformal 3-geometry, 7^ , and initial slicing as specified 
by K. To date, only flat conformal 3-geometries and 
maximal slicings (I will refer to these as the CFMS as- 
sumptions) have been explored extensively. It has been 
pointed out that quasi-equilibrium binary systems in cir- 
cular orbits will not exhibit a conformally flat 3-geometry 
at second post- Newtonian order [j45|. Some consider this 
approximation to be the major source of error in cur- 
rent numerical work to model quasi-equilibrium binary 
configurations and in locating the ISCO0. Certainly, 
the assumption of conformal flatness introduces some er- 
ror. But, it is not clear that this is the most significant 
source of error in these calculations. As stated earlier, I 
suspect that the ad hoc normalization condition used to 
determine VL and construct evolutionary sequences is the 
dominant source of error. It will be important to explore 
the sensitivity of solutions to the choices for 7^ and K 
in order to understand this issue better. 

Of course, it is possible that no regular quasi- 
equilibrium solutions exist when the CFMS assumptions 



are made. As discussed in § IV A, a regular solution 



where f is the usual areal radial coordinate in Kerr-Schild 
coordinates, the spatial metric on the Kerr-Schild slic- 
ing is seen to be conformally flat. Equations (p0|), (p6|), 
and (^6|) are obviously satisfied, and it is easy to verify 
that thequasi-equilibrium boundary conditions given by 
Eqns. ((77]), (|82|), and ( [si]) are also satisfied on the hori- 
zon at f = 2M. Of course, the slicing in this solution 
is not maximal, or even asymptotically maximal. A fea- 
ture of this and other similar non-maximal single hole 
solutions is that a =/= on the horizon. It seems likely 
that this feature will also hold for similar binary configu- 
rations, removing the most obvious obstacle confronting 
the construction of regular quasi-equilibrium solutions. 

It is clear that we will want to explore configurations 
that are not constructed on a maximal slice and we will 
want to consider conformal 3-geometries that are not flat. 
This idea has been explored in similar contexts Jl0|, p"T| ], 
but little has been done other than to demonstrate that 
solutions can be found. The major limitation of these 
works is that the boundary conditions used were not 
well motivated. However, some of the work done in con- 
structing appropriate values for 7^ and K based on the 
Schwarzschild and Kerr geometries should prove useful. 
In the context of quasi-equilibrium solutions, one should 
construct conformal 3-geometries and associated values 



12 



of K that approximate black hole binaries that are at 
rest in the corotating coordinate system. By superpos- 
ing two Schwarzschild holes, one would approximate the 
3-geometry of a pair of corotating black holes. A super- 
position of two Kerr holes could approximate the case 
of irrotational black holes, assuming the spins are cho- 
sen correctly. Again, it will be interesting to asses the 
sensitivity of corotating and irrotational solutions to the 
choice of the conformal 3-geometry. 

An improvement in the choice for 7^ and K would 
come from taking these data from post-Newtonian solu- 
tions for binaries in circular-orbits. An appropriate met- 
ric could be constructed by omitting radiation damping 
terms. One possibility would be to revert the effective 
one-body metric obtained in Ref. |4^| to an appropri- 
ate two body coordinate system, assuming this can be 
done. A metric accurate to 2.5PN order, including spin 
effects, has been given explicitly in Ref. p^ . Another 
possibility is found in Ref. (i^] where a post-Newtonian 
metric is matched to two perturbed Schwarzschild me- 
tircs. However, it is unclear to me whether or not any 



of these solutions are immediately applicable. A usable 
metric must be written in terms of a gauge where the 
lapse and 3-metric smoothly cross the individual black- 
hole horizons. In particular, the lapse should not vanish 
at these horizons. 

Finally, we might consider the use of a slicing and 3- 
metric obtained from a post-Newtonian solution where 
radiation reaction is included. In this case, we could also 
obtain a non-zero solution for u lJ . However, we should 
be cautious in exploring solutions with u y 7^ since this 
violates the principle of quasi-equilibrium. In particu- 
lar, the condition used to determine the orbital angular 
velocity, f2, will most likely not be applicable in this case. 
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